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Abstract 

Boolean networks have been the object of much attention, especially since S. Kauff- 
man proposed them in the 1960's as models for gene regulatory networks. These systems 
are characterized by being defined on a Boolean state space and by simultaneous up- 
dating at discrete time steps. Of particular importance for biological applications are 
networks in which the indegree for each variable is bounded by a fixed constant, as was 
stressed by Kauffman in his original papers. 

An important question is which conditions on the network topology can rule out 
exponentially long periodic orbits in the system. In this paper we consider cooperative 
systems, i.e. systems with positive feedback interconnections among all variables, which 
in a continuous setting guarantees a very stable dynamics. In Part I of this paper we 
presented a construction that shows that for an arbitrary constant < c < 2 and suf- 
ficiently large n there exist n-dimensional Boolean cooperative networks in which both 
the indegree and outdegree of each for each variable is bounded by two (bi-quadratic 
networks) and which nevertheless contain periodic orbits of length at least c". 

In this part, we prove an inverse result showing that for sufficiently large n and for 
< c < 2 sufficiently close to 2, any n-dimensional cooperative, bi-quadratic Boolean 
network with a cycle of length at least c" must have a large proportion of variables 
with indegree 1. Such systems therefore share a structural similarity to the systems 
constructed in Part I. 
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1 Introduction 



In this note we continue our study of the existence of exponentially long periodic orbits in 
bi-quadratic cooperative Boolean systems. The motivation for our interest in this problem 
was described in Part I of this paper [1]. 

For a positive integer n, let [n] := {!,..., n}. An n-dimensional Boolean dynamical 
system or Boolean network is a pair (H, g), where 11 = {0, l}'"] and g : H ^ H. A state 
s{t) at time t will be denoted by s{t) = [si(t), . . . , or simply s = if 

time-dependency is ignored. We will have 

s{t + l)=g{s{t)). (1) 

The cooperative order on 11 is the partial order relation defined by s < r iff Sj < for 
all i £ [n]. The system is cooperative if s{t) < r{t) implies s{t + 1) < r{t + 1). 

We associate a directed graph D with vertex set [n] with the system. A pair < i,j > 
is in the arc set of D iff there exist states s,r £ H such that Si < ri and Sk = r^ for all 
k ji^ i with the property that {g{si))j < {g{ri))j. We will say that the system is quadratic 
if the indegree of all vertices in D is at most two. We call the system b-quadratic if it is 
quadratic and the outdegree of all vertices in D is at most b, where 6 is a positive integer. A 
2-quadratic system is called bi-quadratic. More generally, for positive integers b, r we define 
a (6, r)-Boolean system as a system in which the indegree of all vertices in D is bounded 
by r and the outdegree of all vertices in D is bounded by b. 

In Part I of this paper [3] it was shown that for every < c < 2 and sufficiently large n 
there exist n-dimensional bi-quadratic cooperative Boolean networks that contain periodic 
orbits of length at least c". These systems were constructed by conceptualizing a small 
subset M of the variables as a Turing machine and the set T of the remaining n — |M| 
variables as L circular tapes in such a way that M writes successive codes of the integers 
0, . . . , [c"J on the tape. A judicious choice of coding allowed us to find examples where the 
whole system is cooperative and bi-quadratic. 

While the metaphor of a Turing machine acting on one or several tapes is certainly 
appealing to the human mind, it is an intriguing question whether our construction is, in 
some sense, the only way of producing cooperative bi-quadratic Boolean systems with very 
long cycles. Here we show that this is indeed the case. 

Consider an n-dimensional Boolean system {Ji,g). Then g = [51,... ,gn]- Taking our 
motivation from Boolean models of gene regulatory networks, we call gt the k-th regulatory 
function. If (11,(7) is quadratic, then gi depends on at most two variables ik,jk- If) hi 
addition, {Il,g) is cooperative, and if ik / jk, then we must have g^ = Sj^. A sj^, or gj^ = 
^ik ^ ^jk- I^ these two cases, we will say that gk is strictly quadratic. The only other 
possibility is 5^ = Sj^, ; we will say in this case that gi^. is monic. Note that if g^ is constant 
then we get identical dynamics along attractors if we replace it with the monic function 
9k = •Sfc- Since transient states are irrelevant for our results, we will wlog assume that the 
indegree of each variable is at least one. In a more general setting, we will say that gk is 
monic iff the indegree of k in the digraph D associated with the system is one. 
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Now consider a sequence of variables . . . , km such that Qki^-^ = Sfc^ for ah i £ [m — 1]. 
The dynamics of the system on these variables is analogous to that of a memory tape that 
advances by one position at each time step. A new value may be written to position ki at 
each time step, and this value may be read i time steps later by some regulatory function 
off position /c^+i. If km = ki, the tape is 'read-only,' and a constant regulatory function can 
be considered a special case of a 'read-only' tape of length one. Also, a tape could split into 
two or more branches (say, if 52 = si, 53 = 54 = S2, 55 = S3, qq = S4, etc., then a branching 
would occur at the second variable), but the values on these branches would eventually 
be only copies of each other. Thus any cooperative system that contains monic regulatory 
functions can be conceptualized as a Turing machine acting one or more tapes, possibly 
branching or of varying lengths. This observation motivates the following definition. 

Definition 1 We call an n- dimensional Boolean system an (M, n)-Turing system if at least 
n — M of the regulatory functions are monic. 

While every n-dimensional Boolean system is an (n, n)-Turing system in the sense of the 
above definition, we will use this expression to highlight the fact that if M < n, the roles 
of the 'machine' and the 'tapes' can be neatly separated. Note that we do not require in 
Definition [T] that the system be cooperative. Admittedly, if also monic regulatory functions 
gi^ = -iSjj. may occur in the system, then the connection with the Turing machine metaphor 
becomes more tenuous, but we will still use this terminology as a convenient way to formulate 
our results. 

As indicated above, the systems constructed in [4] are (M(n), n)-Turing systems such 
that lim^^oo = 0. The main result we will prove here shows that if c is sufficiently 

close to 2, then in every bi-quadratic cooperative Boolean system with a periodic orbit of 
length at least c" the 'tapes' must contain the vast majority of all variables. 

Theorem 2 Let a > and let b be a positive integer. Then there exists a positive constant 
c < 2 such that for sufficiently large n, every n-dimensional b-quadratic cooperative Boolean 
system with a periodic orbit of length at least c" is an (an,n)-Turing system. 

We will actually prove a more general result than Theorem [5J The bias A of a Boolean 
function is the fraction of input vectors for which the function outputs 1. Note that the 
strictly quadratic cooperative Boolean functions x A y and x V y have bias A = 0.25 and 
A = 0.75 respectively, whereas monic Boolean functions have bias A = 0.5. There are two 
strictly quadratic Boolean functions with bias A = 0.5, namely the exclusive-or function 
and the equivalence function, but these are not cooperative and cannot occur in cooperative 
Boolean systems. On the other hand, there are cooperative Boolean functions with bias 
A = 0.5 that depend on three input variables; an important example is the function that 
takes the value 1 iff the majority of input variables have the value 1. 

We will say that a Boolean system (11, g) is e-biased if every non-monic regulatory 
function has bias A with |A — 0.5| > e. It follows from the above discussion that cooperative 
quadratic Boolean systems are 0.25-biased, but in general, cooperative Boolean systems 
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with regulatory functions that can take three or more inputs need not be e-biased for any 
e > 0. We will prove the following generalization of Theorem [2j 

Theorem 3 Let e, a > and let b, r he positive integers. Then there exists a positive 
constant c(e, a, b,r) < 2 such that for every c > c(e, a, b, r) and sufficiently large n, every 
n- dimensional e-biased (6, r) -Boolean system with a periodic orbit of length at least is an 
(an,n)-Turing system. 

A Boolean function (p that depends on variables xi, . . . ,xi is canalyzing if there exists 
one input variable Xc, called the canalyzing variable, a Boolean value u called the canalyzing 
value, and a Boolean value v called the canalyzed value such that ip{xi, . . . , xi) = v whenever 
Xc = u. Clearly, all three kinds of regulatory functions that are permitted in quadratic 
cooperative Boolean systems are canalyzing. Another important example of a quadratic 
canalyzing Boolean function is the implication x — > y. It is easy to see that a canalyzing 
Boolean function has bias A = 0.5 iff it is monic. Since there are only finitely many Boolean 
functions on any fixed number of inputs, the following is an immediate consequence of 
Theorem [3l 

Corollary 4 Let a > and let b, r be positive integers. Then there exists a positive constant 
c < 2 such that for sufficiently large n, every n-dimensional (b,r)-Boolean system with a 
periodic orbit of length at least c" , and in which all regulatory functions are canalyzing, is 
an {an,n)-Turing system. 



2 Proof of Theorem [3] 

We will prove Theorem [3] in two stages. In the first stage of the proof we will show that very 
large subsets of the state space 11 of an n-dimensional Boolean system must be balanced in 
a sense that will be defined shortly. In the second stage of the proof we will show that if S is 
the set of states in a periodic orbit of an e-biased (6, r)-Boolean system and S is sufficiently 
balanced, then only a small fraction of the regulatory functions can be non-monic. 



2.1 Balanced subsets of the state space 

Let n = {0, l}!*^' be the state space of an n-dimensional Boolean system. Let S = {s^ : I G 
L} be a sequence of (not necessarily pairwise distinct) elements of 11. If the elements of S 
happen to be pairwise distinct, then we will speak of S being a subset of 11. 
To illustrate the key idea of this section, let i G [n] and consider the ratio 

If 7 > 0, then we will say that S is f3-j-l-balanced if \{i E [n] : \(i{S)—0.5\ > 7}| < /3n. 
More generally, let r S [n] and a : [r] — > {0, 1}. For r-element subsets / = {ii, . . . ,ir} 
of [n] with ii < ■ ■ ■ < ir we define ratios (,f{S) as follows: 
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?/(^) - \L\ 



Define 



CnS) = max{2-l^l - Cf{S) : a G {0, 1}M}. 

If /3, 7 > 0, then we will say that S is P-j-r -balanced if for every family P of pairwise 
disjoint subsets / of [n] with | |J P| > /3n and 1 < |/| < r for each I £ P there exists I £ P 
such that CK*^) < 7- 

We will prove the following. 

Lemma 5 Let r be a positive integer, 7 > and assume 7 < 2"''. Let 

(1 — 2~^ \ / 2"*^ \ ^"'^"'^ 

1-2- + 7J V^^J 

anrf Zei c be a constant such that 

c>2{X{j,r)f. 

Then for sufficiently large n, every subset S of {0, of size > c" is [3-^ -r -balanced. 

Proof: Let /?, 7, r be as in the assumptions, and assume throughout this argument that n 
is a sufficiently large positive integer. Let g > 0, let 1 < c < 2, and let 5 be such that 
1 + Q^<5<\ + 2^7 and 5c^ is an integer. Let us assume that S = {s^ : £ S [5c"]} is 
a sequence of randomly and independently (with replacement) chosen states in {0, of 
length 56^. We will treat and Q as random variables and temporarily suppress their 
dependence on S in our notation. 

Let V £ [r]. For fixed I = {ii, . . . , i^} with ii < ■ ■ ■ < iy and a G {0, l}'"' we define 

a _ 1^1=1 Vie 
- 5c- ' 

where r/J^ = if s^{iu) = o'(u) for all u G [w], and t?^^ = 1 otherwise. 

Clearly, the mean value of r/J is E{r]f) = 1-2-". Note that 2l^l-^J < e iSr]f-E{r]f) > e, 
and hence Cj>eiSrij- E{ri]) > e for at least one a G {0, 1}^ . 

We want to estimate Pr{'qf — E{r]f) > e) for any given fixed e > 0. Note that the 
random variables rjf^ take values in the interval [0, 1] and are independent. This allows us 
to use the following inequality of [8] (see also [12\ [1] for the special case we are considering 
here) . 

Lemma 6 Let Xi,X2, ■ ■ . ,Xm be independent random variables such that < Xj < 1 for 
i G [m] and let X = {Xi + • • • + Xm)/nn. Let n = E[X) and let < e < 1 — fi. Then 
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We will assume until further notice that e < 2^'^ and thus satisfies the assumptions 
of ([2]). Both bounds in ([2]) are of the form A™ for some < A < e~^^ < 1. For the moment, 
assume that A is such such a constant, and let m = 5c^. Now it follows from ^ that 

Pr(7?7 - 1 + 2^"^ > e) < A-^^". 
This implies the following estimate for Q: 

Pr{Ci > e) < 2''A'^'=". 

Now fix < n and consider k pairwise disjoint subsets Ii, . . . ,1^ of cardinality < r each. 
The random variables Q^, . . . , are independent. It follows that 

Pr(Vm G [k] Ci^ >e)< 2^'^A'^^^". 

Let k = (in and let A be the event that there exists family P of pairwise disjoint subsets 
/ of [n] with I U ^1 — P"^ 1 < l-^l < for each I G P such that C/ ^ ^ for each 
I G P. The number of eligible families P is bounded from above by (")^" < n^'^". Thus 
the probability of the event A can be estimated as 

Pr{A) < (2n)'^^"A^"^'=". 
Now note that by Stirling's formula the number of subsets of IT of size c" satisfies 



< — - < -- 

c'- ■ 



Moreover, note that 



Thus for 



r/3 

lim (2n)c" = 1. 

n—^oo 



c > 2A^'^ (3) 

and n sufficiently large, we will have 



(2n)i^A'^^< 



-1 



c / 

This in turn implies that for sufficiently large n and c as in ([3]) 



Pr{A) < (2n)^/5"A/^"^'=" = ({2npX'''Y" < (f£:i)— " < _i_ 



(4) 



Now let us fix £ such that < e < 7. Since 7 < 2 ^, the assumptions of Lemma [6] will 
be satisfied for this choice of e. Let B = B{S) be the set of the first c" pairwise distinct 
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elements of the sequence S, if in fact S has at least c" pairwise distinct elements, and let 
B be undefined otherwise. Let us make a few observations: 

1. Let N = {£ ^ [(5c"] : 31 < j < i = s^} be the number of entries in S that duplicate 
a previous entry. Note that B is defined iff < (6— l)c". In particular, by the choice 
of 6, the set B is defined as long as < g^c^. 

2. Note that the expected value of N can be estimated, for sufficiently large n, fixed 
c < 2, and < g < as 

-— < = o(l)c'^. 

on on ^ ' 



2" 

^e[<5c"] 



In particular, £^(iV) < ^c" 



3. Now it follows from Markov's Inequality 

^-^C" > E{N) > Pr{N >{5- l)c")((5 - l)c'^ 

that for fixed c and sufficiently large n, the set B will be defined with probability 
> 0.5. 

4. Assume B is defined. Observe that for each subset / of [n] and o G {0, Ijtl^ll we have 

&1,,.^S)<&1±1^. (5) 

The first inequality in ([5]) turns into equality if r/^^ = whenever is outside of 5; 
the second inequality in ([5]) turns into equality if r/^^ = 1 whenever is outside of B. 
It follows from the relationship between the 7]J{S)'s and QiS) that 

By choosing g sufficiently close to 0, we can choose 6 as close to one as we need, and 
our choice of e < 7 implies that for 6 sufficiently close to one the inequality Q (S) < e 
will imply the inequality Ci{B) < 7. 

But if there is any subset B of size c" of 11 that is not (3-^-r balanced, then this subset 
will be exactly as likely to be equal to B{S) as any other subset of 11 of the same size. By 
point 3 above, the probability that B{S) exists is greater than 0.5, and thus the probability 
that B{S) exists and is equal to B must be at least 0.5 (^„) . But point 4 above implies 
that if B is not f3-^-r balanced, then B{C) = B implies that the event A has occurred, with 
contradicts inequality ([!]). 
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We derived the contradiction under the assumption that c satisfies inequahty ([3|). Now 
assume 

c>2(A(7,r))^ 

as in the assumption of the lemma. Then we can choose e sufficiently close to 7 and 
A = A(e,r) so that inequality ^ will hold as well for any 6 > 1. By choosing 5 sufficiently 
close to one we will get a contradiction whenever B exists and satisfies > 7- This 

proves Lemma O □ 

2.2 Systems with balanced periodic orbits 

Lemma 7 Let b, r be positive integers, let < e,T < 0.5, let (H, g) be an n- dimensional 
Boolean system, and let S be a periodic orbit of (H, g) . Let k S [n] be such that the bias A 
of gk satisfies |A — 0.5| > e, and let L be the set of input variables of g^- Then either 
C;(5)>^ «K|\}(S)>(l-T)e-5. 

Proof: Assume wlog that A > 0.5 + e; the proof in the case when A < 0.5 — e is symmetric. 
Suppose that < Then there exists a subset 5"* C 5 with \S*\ > (1 — r)|5| such that 
rjJiS*) = 2-1^1 for each a £ {0, 1}^. We conclude that 

\S\Ck = \{seS: Sk = 1}| = \{sGS: gk^s) = 1}\ > \{s G 5* : g^is) = 1}\ 
> (1 - t)|S|A > (1 - t) \S\ (0.5 + e)> \S\ (o.5 + (1 - T)e - ^) , 

and the inequality > (1 — T)e — ^ follows. □ 

Lemma 8 Let {Il,g) be an n- dimensional e -biased {n,r)-Boolean system, letO < t < ^q^j-g, 
let 7 = ^ , 7* = (1 — r)e — ^, and let /3, /?* > 0. Assume S is the set of states in a periodic 
orbit of (Il,g) so that S is both (3 -j-r -balanced and (3* -j* -1-balanced. Then there exists a 
subset J C [n] of size \ J\ < (/5 + r/3*)?i with the property that every non-monic regulatory 
function g^ has at least one input variable in J. 

Proof: Let K = {k £ [n] : C{fc}('S') > 7*}- The assumption on 5 implies that \K\ < f3*n. 

Let Jo be the set of inputs of the variables in K. Then | Jo| < rf3*n. 

Let = [n]\K and let {ki, . . . , kp} C be a set of variables maximal with respect 
to the property that gk^ is non-monic for every q G [p] and the sets Iq of inputs of gk^ are 
pairwise disjoint. Let Ji = Uge[p] ^g- 

By Lemma [7] and the choice of K~^, for each q £ \p] we must have Q (S) > Thus 
the assumption on S implies that | Ji| < l3n. 

On the other hand, by maximality of {ki, . . . ,kp}, every non-monic regulatory func- 
tion gk must have at least one input in the set J := Jq U Ji, and the lemma follows. □ 

Now let {Il,g),£,a,b,r be as in the assumptions of Theorem [3l let 7,7* be as in the 
assumptions of Lemma [51 and assume that > satisfy 
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b 

Let A(7,r), A(7*, 1) be as in Lemma [5l Then we will have 



(6) 



c(e,a,6,r) <max{2(A(7,r))^,2(A(7M)r}. (7) 

To see this, assume n is sufficiently large and S is a periodic orbit of (H, g) of length 
at least c", where c exceeds the right-hand side of ([7]). Then Lemma [5] implies that S is 
/3-7-r-balanced and /3*-7*-l-balanced and thus satisfies the assumptions of Lemma [8l Let J 
be as in the conclusion of Lemma [HI Note that at most b\J\ < an regulatory functions can 
have inputs in J, and it follows that (Jl,g) is an (an, n)-Turing system. 

Note that by the second inequality in ([2|) we will in particular have 

c{e,a,b,r) < max{2e-^^'^ ,2e-^^^'^''^*}. (8) 
This concludes the proof of Theorem [3l □ 

We formulated Theorem [3] as a qualitative result about existence of a constant and 
wrote the proof so as to make the argument as transparent as possible. In the remainder 
of this paper we will use the notation c(e, a, b, r) as shorthand for the largest real number 
for which the conclusion of Theorem [3] holds. 

To arrive at more precise estimates of c(e, a, b, r), we defined 7 = 7* = (1 — T)e — ^, 
and wrote a simple MatLab program for numerically exploring the values of the right-hand 
side of ([7]) for r G (0, j^f^) and (3* E (0, ^). Note that the value of (3 is not a free parameter 
as it is given by Q. 

For the case of bi-quadratic cooperative systems, when e = 0.25 and 6 = r = 2, we 
found an almost perfect linear approximation of the right-hand side of ([7]). 

c(0.25,a,2,2) < max{2(A(7, r))^, 2(A(7* , 1))^*} f« 2 - 0.0041a. (9) 

3 Three additional observations 

Theorem [2] cannot be proved without the assumption of some bound on the outdegrees in 
the associated digraph: note that in the following example the system in question is not 
necessarily bi-quadratic. 

Example 9 Let < c < 2. Then for all sufficiently large n there exists an n- dimensional 
quadratic cooperative Boolean system (H, g) that has only strictly quadratic regulatory func- 
tions and contains a periodic orbit of length c". 

Proof: Let (S, /) be a bi-quadratic cooperative Boolean system of dimension n — 2 that 
contains a periodic orbit of length c"", as constructed in [4]. Let H = {0, let gu = fk 
whenever k < n — 1 and is strictly quadratic, let g^ = si^, A s„ whenever k < n and 
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fk = Sjfc, and let gn-i = dn = Sn-i A s„. Then (H, g') is cooperative, quadratic, and has 
only strictly quadratic regulatory functions. 

Now let s G S be a state in a periodic orbit of length at least c" of (S, /), and define a 
state s* S n by s* = [si, . . . , 5^-2, 1,1]. Then the orbit of s* in (H, g) has the same length 
as the orbit of s in (S, /). □ 

Example 10 Let < c < \/2- Then for all sufficiently large n there exists an n- dimensional 
bi-quadratic cooperative Boolean system (H, g) that has only strictly quadratic regulatory 
functions and contains a periodic orbit of length c^ . 

Proof: If c is as in the assumption, then < 2. Let n be sufficiently large and wlog 
even, and let (S,/) be a bi-quadratic cooperative Boolean system of dimension ^ that 
contains a periodic orbit of length c^", as constructed in [4]. Define an n-dimensional 
Boolean system {Il,g) in such a way that for every monic regulatory function = Sj^, 
we have g^ = fi'n/2+fc = ^ik ^ ^in/2+k ■ This can be done in such a way that (11, g) remains 
bi-quadratic, cooperative, and no new monic regulatory functions are introduced, while the 
dynamics of the system on [^] coincides with the dynamics of (S, /) for all initial states 
with si{0) = s„/2+i(0) for all 

Finally, recall that a subset T = {fci, . . . , k^} of [n] with gt-^^ = Sfc. for all f G [m — 1] 
and gkj = Sfe„ of a Boolean system (Jl,g) can be considered a 'read-only' tape. Let us 
define a generalized read-only tape as a set of variables T = {ki, . . . , k^} ^ [n] such that 
k gT, for all i G [m - 1] we have g^^^, = Sfc, or ^^,+1 = "'Sfc,, and gk^ = Sk^ or 5^1 = "'Sfc™.- 
Let us define the read-only part or strictly monic part R{g) of a Boolean system (11, (7) as 
the union of all its generalized read-only tapes. The following theorem shows that Boolean 
systems with large read-only parts cannot have very long periodic orbits. 

Theorem 11 Let (5 > and let c > 2^~^ . Then for sufficiently large n, no n-dimensional 
Boolean system {Il,g) with \R{g)\ > 5n can have a periodic orbit of length > c". 

Proof: First note that if T = {ki, . . . , k^} is a generalized tape of length m, then the 
dynamics of the system on T is can be described by cyclical shifts, possibly with negations in 
some positions. Thus given any initial state s(0) of the system, the vector [s^.^ (i), . . . , (t)] 
can take at most 2m distinct values throughout the trajectory of s(0). Since R{g) = 
{ii, . . . , in} is the union of pairwise disjoint generalized tapes Ti, . . . , T„ with |Ti| + • • • -|- 
\Ty\ = \R{g)\, it also follows from the same observation that the vector [s^^(t), . . . , s^^(t)] 
can take at most 2/cm ({|ri|, . . . , |r„|}) distinct values throughout the trajectory of s(0), 
where Icm stands for the least common multiple. 

If P{N) denotes the maximum value of Icm ({mi, . . . , mr}) with X]i=i ~ then 
P(N) = ev'^i'^^(i+°(i)) as N ^ 00 (see Chapter 13 of [IT]). It follows that for any 
given initial state s(0), the vector of values of the variables in R{g) can take at most 
2e\/l-^(s')l different values in any periodic orbit. Thus under the assumptions 
of the theorem, the number of different states in any orbit is bounded from above by 
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2(1 <5)n2gV^nTiiri(i+o(i)^ which is less than c" for c as in the assumption and sufficiently 
large n. □ 

4 Discussion 

The dependence of the length of periodic orbits, or limit cycles, in Boolean networks on 
the network architecture has been the subject of several empirical studies. In these studies, 
random networks are drawn from the ensemble of all Boolean networks with prescribed 
restrictions on the architecture, and the dynamics is studied by simulating the trajectories 
for a sample of initial states. Very long limit cycles are a hallmark of the chaotic regime, 
whereas relatively short limit cycles are a hallmark of the ordered regime. 

It has been observed that the dynamics tends to become more ordered and less chaotic 
if the number of inputs for each regulatory functions is small [9] , if the regulatory functions 
are strongly biased [HI [E] , if all regulatory functions are (nested) canalyzing functions [TO] , 
or if there are few negative feedback loops [13] . 

Our approach is different from that of the empirical studies cited above in that we try 
to find conditions on the network architecture that provably exclude periodic orbits whose 
length exceeds certain bounds instead of showing that very long orbits are rare in typical 
networks with certain architectures. Note that bi-quadratic cooperative networks embody 
all the restrictions studied in the cited empirical work: The number of inputs of each 
regulatory function is restricted to two, each regulatory function is (nested) canalyzing, has 
a strong bias unless it is monic, and negative feedback loops are absent. Nevertheless, as 
Part I of this preprint [1] shows, the combination of these properties still does not eliminate 
the theoretical possibility of existence of periodic orbits of length > for any c < 2. 

The novelty of the results presented here is that if, in addition, all or a specified fraction 
of the regulatory functions are required to be non-monic, then there exist constants c < 2 
such that the resulting Boolean networks cannot have any periodic orbits of length > c"' 
whatsoever. As Theorem [3] shows, this will be true as long as there are fixed bounds on 
the number of inputs and outputs of each regulatory function as well as on their biases, 
regardless of the prevalence of negative feedback loops or whether the regulatory functions 
are actually canalyzing. 

One can interpret Theorem [2] in a different way. As noted earlier, variables with monic 
regulatory functions in cooperative Boolean systems just record the values of other variables 
at some time in the past (less than n steps earlier). Thus if we allow time delays in the 
definitions of regulatory functions, we can remove all but the first variable on each 'tape' 
and define a Boolean delay system on the remaining variables that will have equivalent 
dynamics, in particular, that will have periodic orbits of the same length as the original 
system. 

Formally, let d be a positive integer and consider the state space 11 of an M-dimensional 
Boolean system and a map h so that h : U^'^ 11. Then (11, h) defines an M-dimensional 
Boolean delay system with maximum delay d whose dynamics is given, for all times t > d 
by 
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s{t + 1) = h{s{t),s{t - 1), s{t-d+ 1)). 



(10) 



We associate to /i a function H : 11^^ UJ-'^ defined by H{s) = [s'^ , . . . , s'^ , h{s)] for 
s = [s^ , . . . , s'^]. We will say that an n-dimensional Boolean system (T,,g) is induced by 
an M -dimensional Boolean delay system (IT, h) with maximum delay d if there exists an 
injection F = (F1,F2) : [n] [M] x [d] such that for the map /* : nl"^! S defined by 

{f*{^))i = ^Fifi) ^'^y f '■ ^ ^ with f*of = idwe have 



In particular, as we indicated above, a cooperative n-dimensional Boolean system 
will be induced by a Boolean delay system with maximum delay n whose variables are 
only the ones with non-monic regulatory functions /fc, together with one variable from each 
'read-only' tape. This allows us to reformulate Theorem [2] as follows: 

Theorem 12 Let a > and let b be a positive integer. Then there exists a positive constant 
c < 2 such that for sufficiently large n, every n-dimensional b-quadratic cooperative Boolean 
system with a periodic orbit of length at least c"" is induced by a Boolean delay system with 
maximum delay at most n and dimension at most an. 

It follows that all cooperative bi-quadratic Boolean systems with very long periodic 
orbits must be induced by Boolean delay systems with the same properties but of much 
smaller dimension. 

Boolean delay systems as defined above are a special case of continuous-time Boolean 
delay systems as introduced in [21 [5] . See [6] for a comprehensive survey and additional 
references. They may be relevant as a modeling tool for gene regulatory networks, since 
gene regulation always involves a delay between gene transcription and the time when the 
translated gene product becomes available as a regulator, such as a transcription factor. 
In fact, they are closely related to the conceptual framework of 'kinetic logic' that was 
developed by R. Thomas in [141 [TS] specifically for the study of gene regulation. 

Our results appear to us of some philosophical interest given the metaphors discussed 
here and in the introduction. Our theorems have the advantage over numerical results of 
giving bounds that are rigorous and universally valid. We want to emphasize though that 
our proof of Theorem [3] only shows existence and gives an upper bound for c(e, q,6, r); we 
do not know whether this estimate is anywhere close to optimal. 

For example, 



c(0.25,0.1,2,2) < 1.9996. 

Note that if a bi-quadratic cooperative system of sufficiently large dimension n has a 
periodic orbit of length > c(0.25, 0.1, 2, 2)" , then at least 90% of all regulatory functions 





c(0.25, 1,2,2) < 1.9959, 
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must be monic; if such a system has a periodic orbit of length > c(0.25, 1,2,2)"', then at 
least some of the regulatory functions must be monic. On the other hand. Example [10] gives 
a lower bound 

v/2 < c(0.25, 1,2,2). 

It will be an interesting direction for future research to narrow the gap between the 
upper and lower bounds. 
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